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Abstract: We use data from Type la Supernovae (SNIa), Baryon Acoustic Oscillations 
(BAO), and Cosmic Microwave Background (CMB) observations to constrain the recently 
proposed teleparallel dark energy scenario based on the teleparallel equivalent of General 
Relativity, in which one adds a canonical scalar field, allowing also for a nonminimal cou- 
pling with gravity. Using the power-law, the exponential and the inverse hyperbolic cosine 
potential ansatzes, we show that the scenario is compatible with observations. In partic- 
ular, the data favor a nonminimal coupling, and although the scalar field is canonical the 
model can describe both the quintessence and phantom regimes. 
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1 Introduction. 

Recent cosmological observations [1-5] indicate that the observable universe experiences 
an accelerated expansion. Although the simplest way is the consideration of a cosmological 
constant, there are two alternative ways to explain this behavior. The first is to modify 
the content of the universe by introducing the dark energy sector, which can be based on a 
canonical scalar field (quintessence) [6-22], a phantom field [23-28], or on the combination 
of both these fields in a unified scenario called quintom [29-35] (see [36] for a review). 
The other approach is to modify the gravitational sector itself (see [37] for a review and 
references therein). 

The teleparallel dark energy is a recently proposed scenario that tries to incorporate 
the dark energy sector [38] . It is based on the "teleparallel" equivalent of General Relativity 
(TEGR), that is on its torsion instead of curvature formulation [39, 40], in which one adds 
a canonical scalar field, and the dark energy sector is attributed to this field. In the case 
where the field is minimally coupled to gravity the scenario is completely equivalent to the 
standard quintessence [6-11], both at the background and perturbation levels. However, in 
the case where one switches on the nonminimal coupling between the field and the torsion 
scalar, that is the only suitable gravitational scalar in TEGR, the resulting scenario has 
a richer structure, exhibiting quintessence-like or phantom-like behavior, or experiencing 
the phantom-divide crossing [38]. The richer structure is manifested in the absence of a 
conformal transformation to an equivalent minimally-coupled model with transformed field 
and potential, which is known to be able to describe only the quintessence regime. 

Since the teleparallel dark energy exhibits interesting cosmological behavior, in the 
present work we use observations in order to constrain the parameters of the model. In 
particular, we use data from Type la Supernovae (SNIa), Baryon Acoustic Oscillations 
(BAO), and Cosmic Microwave Background (CMB) to plot likelihood-contours for the 
present dark-energy equation-of-state, matter density parameter and nonminimal coupling 
parameters, respectively. 
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The paper is organized as follows: In Section 2 we present the scenario of the teleparallel 
dark energy and derive the relevant equations for the cosmological evolution. In Section 3 
we use observational data to produce likelihood-contours of the model parameters. Finally, 
Section 4 is devoted to the summary of our results. 



2 Teleparallel Dark Energy 

Let us briefly review the teleparallel dark energy. As we stated in Introduction, it is based 
on the "teleparallel" equivalent of General Relativity (TEGR) [39, 40], in which instead 
of using the torsionless Levi-Civita connection one uses the curvatureless Weitzenbock 
one. The dynamical objects are the four linearly independent vierbeins (these are parallel 
vector fields, referred to as the appellations "teleparallel" or "absolute parallelism"). It 
is interesting to note that the torsion tensor is formed solely from products of the first 
derivatives of the tetrad. 

In particular, the vierbein field ga{x^) forms an orthonormal basis for the tangent space 
at each point x M , that is • e# = t]ab, where t\ab = diag(l, —1,-1, —1), and furthermore 
the vector can be analyzed with the use of its components in a coordinate basis, 
that is eA = e^d^ 1 . In such a construction, the metric tensor is obtained from the dual 
vierbein as 

9nu{x) = f]AB e^(x) (x). (2.1) 

W A 

Consequently, the torsion tensor of the Weitzenbock connection T v/1 [43] reads 

T%, = - = e\ (d^ - d u e£), (2.2) 

W A 

where F v/Jt = e\ <9 M e~. 

In the present formalism, all the information concerning the gravitational field is in- 
cluded in the torsion tensor T^ u . As described in [40], the "teleparallel Lagrangian" can 
be constructed from this torsion tensor under the assumptions of invariance under gen- 
eral coordinate transformations, global Lorentz transformations, and the parity operation, 
along with requiring the Lagrangian density to be the second order in the torsion tensor. 
In particular, it is the torsion scalar T, given by [39, 40, 44, 45]: 

C = T = -^ v T pilv + \t^ v T vw - T pl fT v %. (2.3) 
The simplest action in a universe governed by teleparallel gravity is 

(2.4) 

where e = det(e^) = y/—g (one could also include a cosmological constant). Variation with 
respect to the vierbein fields provides equations of motion, which are exactly the same as 





' T 


j d A xe 


-z~n + *-"m 
2k z 



1 We follow the notation of [41, 42], that is Greek indices n,is,... and capital Latin indices A,B,... run 
over all coordinate and tangent space-time 0, 1, 2, 3, while lower case Latin indices (from the middle of the 
alphabet) i,j, ... and lower case Latin indices (from the beginning of the alphabet) a,b,... run over spatial 
and tangent space coordinates 1, 2, 3, respectively. Finally, we use the metric signature (+, — , — , — ). 
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those of GR for every geometry choice, and that is why the theory is called "teleparallel 
equivalent to General Relativity" . 

In principle one has two ways of generalizing the action (2.4), inspired by the corre- 
sponding procedures of the standard General Relativity. The first is to replace T by an 
arbitrary function f(T) [41, 42, 46-72], similar to f(R) extensions of GR, and obtain new 
interesting terms in the field equations. The other, on which we focus in the present work, 
is to add a canonical scalar field in (2.4), similar to the GR quintessence, allowing for a 
nonminimal coupling between it and gravity. This field will constitute the dark energy 
sector, and thus the corresponding scenario is called "teleparallel dark energy" [38]. In 
particular, the action will simply read: 



S 



d xe 



V{</>) + £ r , 



(2.5) 



Concerning the nonminimal coupling we emphasize that the nonminimal coupling will be 
between the torsion and the scalar field (similar to the standard nonminimal quintessence 
where the scalar field couples to the Ricci scalar). 

Variation of the action (2.5) with respect to the vierbein fields yields the equations of 
motion 

2 + 2^ 2 ) e-%(ee? A S p n - e A T% x S p ^ - \e\T 



- v(<t>) 



+ e fi A d l/ <f>d„<f> 



+A& A Sr<t>(d^) = e A T p », 



(2.6) 



where T p v stands for the usual energy-momentum tensor. Therefore, for a flat Friedmann- 
Robertson- Walker (FRW) background metric 



ds 2 = dt 2 - a 2 (t) 5ijdx l dx j 



(2.7) 



and a vierbein choice of the form 



-A 



e',; = diag(l, a, a, a) 



(2.8) 



where t is the cosmic time, x l are the comoving spatial coordinates and a(t) is the scale 
factor, we obtain the corresponding Friedmann equations: 



H 2 
H 



n 2 ( 
TV 



K 

"2~ 



+ Pr. 



+ P<j> + Pm+Pn 



(2.9) 
(2.10) 



where H = a/ a is the Hubble parameter and a dot denotes differentiation with respect to t. 
In these expressions, p m and p m are the matter energy density and pressure, respectively, 
following the standard evolution equation p m + 3H(l + w m )p m = 0, with w m = p m /Pm the 
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matter equation-of-state parameter. Additionally, we have introduced the energy density 
and pressure of the scalar field 

P<t> = \4> 2 + v{4>)-^H 2 4> 2 , (2.11) 

V4> = \j> 2 - V{4>) + ^H4>4> + £ (3H 2 + 2Hj </> 2 . (2.12) 

Moreover, variation of the action with respect to the scalar field provides its evolution 
equation, namely: 

4> + 3H<j) + 6£H 2 (l) + V'(<l)) = 0. (2.13) 

Note that in the above expressions we have used the useful relation T = —6H 2 , which 
straightforwardly arises from the calculation of (2.3) for the flat FRW geometry. 

In this scenario, similar to the standard quintessence, dark energy is attributed to the 
scalar field, and thus its equation-of-state parameter (wde) reads: 

w DE = w^ = ^. (2.14) 

P<t> 

As a result, one can see that the scalar field evolution (2.13) leads to the standard relation 

^ + 3^(1 + ^)^ = 0. (2.15) 

The teleparallel dark energy proves to exhibit a very interesting cosmological implica- 
tion [38, 73]. In the minimally-coupled case the cosmological equations coincide with those 
of the standard quintessence, both at the background and perturbation levels. However, 
when the nonminimal coupling is switched on, one can obtain a dark energy sector being 
quintessence-like, phantom-like, or experiencing the phantom-divide crossing during evo- 
lution, a behavior that is much richer comparing to General Relativity (GR) with a scalar 
field [38]. Therefore, it is both interesting and necessary to use observations to constrain 
the parameters of the scenario. This is performed in the next section. 

3 Observational Constraints 

We use Type la Supernovae (SNIa) from the Supernova Cosmology Project (SCP) Union2 
compilation [74], Baryon Acoustic Oscillations (BAO) data from the Two-Degree Field 
Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey Data Release 7 (SDSS 
DR7) [75], and the Cosmic Microwave Background (CMB) radiation data from Seven- 
Year Wilkinson Microwave Anisotropy Probe (WMAP) observations [76] to examine the 
teleparallel dark energy scenario. Since the model includes the scalar-field potential, we 
perform our analysis for three different potential cases, given by: 

• Power-Law potential 

This potential class is common in cosmology [77-81]. Although we can straightfor- 
wardly perform our analysis for an arbitrary exponent, for simplicity we focus on the 
most interesting quartic case 

V(4>) = (3.1) 
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• Exponential potential 

Exponential potentials are also very common in the literature [82-85] , which are nec- 
essary to be considered in every observational constraining analysis. In the following 
we use the ansatz 

V(<f>) = V e~ KX *. (3.2) 

• Inverse hyperbolic cosine potential 

The use of hyperbolic cosine potential, or power-law functions of it, has also many 
cosmological implications [86-88]. Although we could perform our analysis for an 
arbitrary exponent, in the following we focus on the inverse case, namely: 

cosh(K<pj 

We examine the constraints on the model parameters and the present values of the 
density parameters, following the x 2 -method for the recent observational data. The detailed 
analysis method for SNIa, BAO and CMB data is summarized in Appendix. In general, we 
are interested in producing the likelihood contours for physically- interesting parameters, 
namely the present dark-energy equation-of-state parameter wde i the present matter 
density parameter £l m o and the nonminimal coupling parameter £. We mention here that £ 
must always be bounded according to a physical constraint, namely it must lead to positive 
Pde and H 2 in relations (2.9) and (2.11), respectively. In practice, £ is found to be mainly 
negative (in our convention), and only a small window of positive values is theoretically 
allowed. In our analysis, for each of the three potentials, we fit three parameters, namely 
wde , the dimensionless Hubble parameter h and fimo(£)> an d then we draw the likelihood- 
contours for Id and 2a confidence levels. 

In Fig. 1 we present the likelihood contours for wde and fl m o with the teleparallel 
dark energy scenario under the quartic potential (3.1). As we observe, the scenario at hand 
is in agreement with observations, and as expected, it can describe both the quintessence 
and phantom regimes. Since the scalar field is canonical, it is a great advantage of the 
present model. 

In Fig. 2 we present the likelihood contours for wde and the nonminimal coupling 
parameter £, for the quartic potential (3.1). Interestingly enough we observe that the 
nonminimal coupling is favored by the data, and in particular a small £ is responsible for 
the quintessence regime, while a larger one leads to the phantom regime. Note that the 
best-fit values of WDE \b.f ~ —0.98 is very close to the cosmological constant. 

In Fig. 3 we present the likelihood contours for wde and fJ m o, for the teleparallel 
dark energy scenario under the exponential potential (3.2). As we observe, this scenario 
is consistent with observations, and it can describe both the quintessence and phantom 
regimes, with the phantom regime favored by the data. Furthermore, in Fig. 4 we present 
the likelihood contours for wde and £, for the exponential potential (3.2). From this graph 
we deduce that a non-minimal coupling is favored by the data, and we observe that wde - 
values close to the cosmological constant bound, either above or below it, can be induced 
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Figure 1. Contour plots of the present dark-energy equation- of- state parameter wde q versus 
the present matter density parameter Q m o under SNIa, BAO and CMB observational data in the 
teleparallel dark energy scenario with the quartic potential V{4>) = Vb</> 4 . The curves correspond to 
la and 2a confidence levels, respectively, and the cross marks the best-fit point. 

Wdeo 



-0.6 
-0.8 
-1 
-1.2 



-0.6 -0.5 -0.4 -0.3 

% 

Figure 2. Contour plots of the present dark-energy equation- of- state parameter wde versus the 
nonminimal coupling parameter £ under SNIa, BAO and CMB observational data, in the teleparallel 
dark energy scenario with the quartic potential V((f>) = Vo0 4 . The curves correspond to la and 2a 
confidence levels, respectively, and the cross marks the best-fit point. 

by a relative large ^-interval, which is an advantage of this scenario. It is interesting to 
mention that the exponential potential was used as an explicit example in [38], and our 
current analysis verifies its results. 





-6- 



Wdeo 
-0.9 i 




Qlllo 



Figure 3. Legend is the same as Fig. 1 but with the exponential potential V(<fi) = V$e kA( ^. 
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Figure 4. Legend is the same as Fig. 2 but with the exponential potential V((f>) = Voe~ KX ^ . 

In Fig. 5 we depict the likelihood contours for wde and fi m o> under the inverse hyper- 
bolic cosine potential (3.3). As we can see, this scenario is consistent with observations. 
However, if we desire to avoid divergences in the wde evolving history, we are restricted 
in the phantom regime. In addition, the best-fit value of WDE \b.f ~ —1.02 is very close to 
the cosmological constant. Furthermore, in Fig. 6 we present the corresponding likelihood 
contours for wde and £. Similarly to the previous cases we can see that the nonminimal 
coupling is favored by the data. It is interesting to note that wde -values close to the 
cosmological constant bound can be induced by a relative large ^-interval. 
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Figure 5. Legend is the same as Fig. 1 but with the inverse hyperbolic cosine potential V(<fi) = 
Vq I cosh(n4>) . 
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Figure 6. Legend is the same as Fig. 2 but with the inverse hyperbolic cosine potential V{4>) = 
Vq j cosh(n(f>) . 

We close this section with a comment on the positive values of the nonminimal cou- 
pling £. As we mentioned above, the positivity requirement for pde and H 2 leads £ to be 
negative, with only a small window of positive values theoretically allowed. Now, in prac- 
tice, if we perform our fitting procedure in the positive £ region for the inverse hyperbolic 
cosine potential, we find that positive £ is excluded. However, for the quartic and expo- 
nential potentials we find the interesting result that for the theoretically allowed positive 
£ (0 < £ < 0.2) wde is always very close to a constant wde with \wde ~ Wde\ ^5 1CT 3 . 
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The reason is that the scenario of the teleparallel dark energy for positive £ (sufficiently 
small in order for the positivity of pde arid H 2 not to be spoiled) always results in the 
stabilization of u>de close to the cosmological constant value, as can be proven by a de- 
tailed phase-space analysis [89]. Such a behavior is an advantage from both observational 
and theoretical point of view. We would like to note for comparison that, in the case of 
quintessence in the Einstein gravity and in the presence of a dust-like component (CDM 
and baryons), dark energy with an exactly constant wde > — 1 is possible for the inverse 
hyperbolic sine potential in some power as was first independently shown in Refs. [90] 
and [91]. 

4 Conclusions 

In the present work we have used observational data to impose constraints on the pa- 
rameters of the teleparallel dark energy scenario [38], which is based on the teleparallel 
equivalent of General Relativity (TEGR), that is on its torsion instead of curvature formu- 
lation [39, 40]. In this model one adds a canonical scalar field, in which the dark energy 
sector is attributed, allowing also for a nonminimal coupling between the field and the 
torsion scalar. Thus, although the minimal case is completely equivalent with the standard 
quintessence, the nonminimal scenario has a richer structure, exhibiting the quintessence- 
like or phantom-like behavior, or experiencing the phantom-divide crossing [38]. 

In particular, we have fitted data from Type la Supernovae (SNIa), Baryon Acoustic 
Oscillations (BAO), and Cosmic Microwave Background (CMB) observations to constrain 
the present dark-energy equation-of-state parameter wde , the present matter density 
parameter fi m o and the nonminimal coupling parameter £. Furthermore, in order to be 
general, for the scalar-field potential we have taken three ansatzes, namely the power-law, 
the exponential and the inverse hyperbolic cosine ones. 

For the power-law (quartic) potential we have seen that teleparallel dark energy is 
compatible with observations and, as expected, it can describe both the quintessence and 
phantom regimes. Additionally, we have shown that the negative nonminimal coupling is 
favored by the data, and in particular a small £ is responsible for the quintessence regime, 
while a larger one leads to the phantom regime. For the exponential potential we have 
demonstrated that both the quintessence and phantom regimes can be described, with the 
phantom regime favored by the data. Moreover, we have found that a negative non-minimal 
coupling is favored and we have observed that wde -values close to the cosmological con- 
stant bound, either above or below it, can be induced by a relative large ^-interval, which 
is an advantage of this scenario. For the inverse hyperbolic cosine potential we have shown 
that wde is restricted in the phantom regime, while £ is restricted to negative values, with 
a relatively large £- interval being able to lead to wde -values close to —1. We remark that 
positive values of £ are excluded for the inverse hyperbolic cosine potential, while for the 
power-law and exponential ones wde is very close to —1. 

In summary, the scenario of the teleparallel dark energy is compatible with observa- 
tions, for all the examined scalar-field potentials. Furthermore, although the scalar field 
is canonical, the model can describe both the quintessence and phantom regimes. These 
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features are an advantage from both observational and theoretical point of view, and they 
make the scenario at hand a good candidate for the description of nature. Finally, the data 
favor a nonminimal coupling, and thus the model is distinguishable from the standard 
quintessence, since the two scenarios are equivalent only for the minimal coupling. 

An interesting and necessary investigation would be to go beyond the background 
analysis, and examine observables that arise at the perturbation level, such are the growth 
of matter overdensities and the gravitational-wave spectrum, since these could also clarify 
possible Lorentz- violation problems that are not seen at the background level [92-94] . Since 
this study lies beyond the scope of the present manuscript, it is left for future investigation. 
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A Analysis method for the observational data 

In this appendix, we explain the methods for the elaboration of observational data from 
Type la Supernovae (SNIa), Baryon Acoustic Oscillations (BAO) and Cosmic Microwave 
Background (CMB) radiation. The % 2 of the combined observational data is given by 

X 2 = Xsn + xIao + Xcmb- (A- 1 ) 

In our fitting procedure we use the simple x 2 method, rather than the Markov-chain Monte 
Carlo (MCMC) approach such as CosmoMC [95]. In the following we describe the calcula- 
tion for the various xf °f each observational dataset (for detailed explanations on the data 
analysis see e.g. [96, 97]). 

a. Type la Supernovae (SNe la) 

SNe la observations provide the information on the luminosity distance Dl as a func- 
tion of the redshift z. The theoretical distance modulus fj,th is defined by 

IHh(zi) = 51og 10 D L (zi) + hq, 

where fi = 42.38 - 51og 10 /t, with h = iJ /100/[km sec" 1 Mpc" 1 ] [76]. The Hubble-free 
luminosity distance for the flat universe is described as 

. . , . r* dz' 

where E(z) = H(z)/Hq, with 

W. = yft® (1 + zf + Oj 0) (1 + z f + fiW (1 + ~fM^*) . 



- 10 - 



Here, = (1 + 0.2271-/V e ff ), where is the present fractional photon energy 
density and iV e g = 3.04 is the effective number of neutrino species [76]. We mention that 
H(z) is evaluated by using numerical solutions of the Friedmann equation. 
The x 2 °f the SNe la data is given by 

2 _ IP-abaiZi) ~ fJ-th(zi)] 2 / a n 

XSN — /-^i ^2 ' \ A - Z ) 

i % 

where /z Q bs is the observed value of the distance modulus. In the following, subscriptions 
"th" and "obs" denote the theoretically predicted and observed values, respectively. xIn ^ s 
minimized with respect to /xo, which relates to the absolute magnitude, since the absolute 
magnitude of SNe la is not known. xIn m (A-2) is expanded as [98] 

,2 



with 



Xsn = A - 2fi B + (4C , 



_ [Mobs(gi) — Hth(Zi', Ho 



Thus, the minimum of x| N with respect to fiQ is expressed as 

Xsn = A - ^ . (A.3) 

In our analysis we apply expression (A.3) for the x 2 minimization and we use the Super- 
nova Cosmology Project (SCP) Union2 compilation, which contains 557 supernovae [74], 
ranging from z = 0.015 to z = 1.4. 

b. Baryon Acoustic Oscillations (BAO) 

The distance ratio of d z = r s (za)/Dy(z) is measured by BAO observations, where Dy 
is the volume-averaged distance, r s is the comoving sound horizon and z& is the redshift at 
the drag epoch [75]. The volume- averaged distance Dy(z) is defined as [99] 



D v (z) 



(l + z) 2 D\{z) 



1/3 



H(z). 

where Da(z) is the proper angular diameter distance for the flat universe, defined by 

1 f z dz' 



H(z>) 
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The comoving sound horizon r s (z) is given by 

I /•!/(!+») 

r s (z) 



VZJo 



da 



a 2 H{z' = l/o-l)*/l+ (3^ 0) /4^ 0) 



where 



(0) 



2.2765 x 10- 2 h~ 2 and n. 



(0) 



2.469 x 10 5 h 2 are the current values of 



baryon and photon density parameters, respectively [76]. The fitting formula for z& is 
given by [100] 



2d 



i29i(n^h 2 ) - 251 

1 + 0.659(fiSP /i 2 ) - 828 



i + M^V 



62 



with 



o.3i3(o^ 2 )-°- 419 

0.238 «/i 2 ) ' 223 



1 + 0.607 (O^ 2 ) - 674 



The typical value of z^ is about 1021 for fifj' = 0.276 and h = 0.705. 

According to the BAO data from the Two-Degree Field Galaxy Redshift Survey (2dF- 
GRS) and the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) [75], the distance 
ratio d z at two redshifts z = 0.2 and z = 0.35 is measured to be d°=o 2 = 0.1905 ± 0.0061 



and ^0.35 



0.1097 ± 0.0036, with the inverse covariance matrix: 



C 



1 

BAO 



30124 -17227 
-17227 86977 



Finally, the x 2 for the BAO data is calculated as 



Xbao 



„th 

^,BAO 



-obs 
c i,BAO 



BAO, 



„th 

j,BAO 



X 



X 



„obs 

j,BAO ) > 



where x i>B AO = (do.2,d . 



35 J 



c. Cosmic Microwave Background (CMB) radiation 



The CMB observational data are sensitive to the distance to the decoupling epoch 
z* [76] . Hence, by using these data we obtain constraints on the model in the high redshift 
regime [z ~ 1000). 

The acoustic scale I a and the shift parameter 1Z [101] are defined by 



Ia(z*) = (1 + z*) 



7tDa(z*) 
r s (z*) 

ft(z*) = ^Q^H (1 + z,) £>a(«0, 
where z* is the redshift of the decoupling epoch, given by [102] 



= 1048 



1 + 0.00124 [n ( ° ] h 2 



-0.738 



i + gi (f4°v) 



,'/2 
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with 

/ , n , \ -0.238 

0.0783 (nrh 2 



/b n J _ 0.560 

■91 - x 0.763' 92 - — — TTIT- 

l + 39.5(f>[ 0) /i 2 J 1 + 21.1« 

We use the data from Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) ob- 
servations [76] on CMB. 

The x 2 of the CMB data is 

XCMB = ( x i|cMB ~~ x ?,CMb) (Comb) ij ( x j!cMB ~~ x ffiMB 

where x^cmb = {Ia(z*),'R'(z*)i z*) an d Cqmb ^ s ^ e inverse covariance matrix. The data 
from WMAP7 observations [76] lead to l A (z*) = 302.09, H{z*) = 1.725 and = 1091.3 
with the inverse covariance matrix: 

/ 2.305 29.698 -1.333 
29.698 6825.27 -113.180 
-1.333 -113.180 3.414 



n-i 

°CMB 



V 
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